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ABSTRACT 


Six space vehicle configurations are investigated to determine 
their relative merits as pertaining to rigid body control characteristics 
and structural bending frequencies. Rigid body theory was used to 
illustrate the factors which govern the stability of a vehicle in flight. 

The aerodynamic moment coefficient (C^), which is one of these factors, 
was determined from slender body theory. The thrust moment co- 
efficient (C 2 ) was determined for each vehicle and the ratio C 2 /C 1 was 
chosen as the control- stability parameter to be used in the compar ison 
of the vehicles. The bending frequency parameter chosen was 
It was observed that long slender vehicles exhibit better control charac- 
teristics than short vehicles with small length-to-diameter ratios, but 
the short vehicles, in turn, have higher bending frequencies. For conical, 
single-diameter cylindrical, and dual-diameter cylindrical vehicles with 
the same approximate length, the conical vehicle exhibits better bending 
and control characteristics than the other two, while the single-diameter 
cylindrical is better than the dual-diameter cylindrical vehicle. 

An arbitrary equation was developed for the weighted diameter, 

D w , of the vehicles, and a vehicle shape parameter was chosen. This 
parameter varies nearly linearly with the structural bending frequencies 
of the vehicles. 
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radius of propellant, 
height of propellant in tank, 
gain factor of attitude channel, 
coefficient of control damping, 
thrust of engine. 

effective moment of inertia of total vehicle about its C. G. 
radius of gyration of vehicle about C. G. 
mass of vehicle, 
coordinate of swivel point. 

coordinate of center of instantaneous rotation, 
bending frequency of the nth mode, 
aerodynamic restoring moment coefficient, 
thrust moment coefficient, 
center of pressure, 
angle of attack at C. G. 

deflection of swivel engine relative to center line of vehicle. 

dynamic pressure. 

diameter of vehicle airframe. 

diameter of base of vehicle 

D/D 

o 

normal lift coefficient, 
length of vehicle. 

average thickness of skin of vehicle, 
weighted diameter of vehicle, 
length of nth section of vehicle, 
stiffness of vehicle, 
mass distribution. 
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by Nathan L. Beard 
SUMMARY 


A vibration analysis of six space vehicle configurations was 
conducted to determine the first three bending modes and their cor- 
responding frequencies. Rigid body theory was used to illustrate the 
factors which govern the stability of a vehicle in flight. The aerodynamic 
moment coefficient (C|), which is one of these factors, was determined 
from slender body theory. The thrust moment coefficient (C 2 ) was 
determined for each vehicle and the ratio C 2 /Cj was chosen as the 
control-stability parameter to be used in the comp ariso n of the vehicles. 
The bending frequency parameter chosen was f n /VfC]J . 

It was found that vehicles with large L/D ratios exhibited higher 
C 2 /C 1 values, or better control characteristics than vehicles with 
small L/D ratios; but the short vehicles, in turn, have higher bending 
frequencies. For conical, single-diameter cylindrical, and dual- 
diameter cylindrical vehicles with the same approximate length, the 
conical vehicle exhibits more favorable bending and control char- 
acteristics than the other two, while the single-diameter cylindrical 
is better than the dual-diameter cylindrical vehicle. 

An arbitrary equation was developed for the weighted diameter, 

D w , of the vehicles, and a vehicle shape parameter was chosen. This 
parameter varies nearly linearly with the structural bending frequencies 
of the vehicles. 
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I. INTRODUCTION 


This investigation was conducted primarily for the purpose o'f 
obtaining information relative to the configur ative design of space 
vehicles with emphasis placed on obtaining larger structural bending 
frequencies. It is desirable for space vehicles to have large natural 
bending frequencies from the stability point of view. For instance, 
when the first bending mode frequency is approximately the same as 
the control or fuel sloshing frequencies, interaction between their 
vibrational modes can bring about control saturation and possibly the 
destruction of the vehicle; therefore, it becomes necessary to vary 
certain parameters to provide a frequency gap between these modes. 

The control frequency may be varied by the slection of different gain 
factors a Q , a^, b Q , and g£ in the control circuit, which are functions 
of time only, while the sloshing frequency can be varied by the proper 
choice of the tank form and the height -to-radius ratio, h/a, where h 
is the height of the fuel in the tanks. With these considerations in 
mind, it is evident that the space vehicle with the higher first mode 
natural bending frequency offers a more expanded or wider range from 
which to choose the gain factors and h/a parameter. 

The stability of a vehicle from the aerodynamic point of view 
depends upon the magnitude and sign of C^/Cj. The more positive this 
ratio the better is the stability of the vehicle. The thrust moment 
coefficient, generally, is positive; however, when the center of pressure 
is forward of the center of gravity of the vehicle, is negative. This 
indicates that a destabilizing moment acts about the center of gravity of 
the vehicle. 

The following analysis of six different vehicle configurations 
shows their favorable and unfavorable qualities pertaining to the 
parameters mentioned above. 
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II. DESCRIPTION 


Six vehicle configurations were considered in this analysis. Five 
of these were chosen as reasonable designs for space vehicles, while 
the Saturn C-l was selected as the standard in this study. Configura- 
tions 2, 3, and 5 were chosen with the same base diameter and approx- 
imate length; however, the overall shape of each was varied to determine 
the effects this would have on the aerodynamic stability and the structural 
bending frequencies. One large L/D and one small L/D vehicle was 
chosen {Vehicles 1 and 4) to contrast the stability and bending char- 
acteristics. The Saturn C-l was chosen, as configuration 0, to be 
compared with the other five since it had at least one characteristic 
of each vehicle. Each of the vehicles was idealized as follows: 

1. The vehicles empty weight was derived only from the 
structural weight of its skin. 

2. The skin thickness of each*vehicle was assumed to be 
constant throughout its length. 

3. All of the vehicles were assumed to be constructed from 
an aluminum alloy with a density of 0. 0975 pounds per 
cubic inch. 

4. The length, diameter, and thickness of the five vehicles 
chosen as reasonable designs were determined so as to 
provide a structural weight and internal capacity equal to 
that of the Saturn vehicle C-l. 

5. The nose cones were given a slope identical to that of the 
Saturn Vehicle (15°, half vertex angle). 

6. Each vehicle was filled completely with water to simulate 
fuel. 

Hereafter the Saturn vehicle will be referred to as model number 
0 while the other five configurations will be referred to as model numbers 
1, 2, 3, 4, and 5. These configurations are shown in Figure 2. 

Since the internal volume and structural weight of each vehicle 
was to be identical to that of vehicle number 0, their length and thickness 
could be easily calculated from the equations for the volume and surface 
areas of a cone and cylinder. This data is given in Table I. 
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Next, the mass distributions, empty and full were plotted as 
shown in Figures 11 through 16. By using a modified Stodola method 
(c. f. Reference 3) in conjunction with the IBM 7090 computer, the 
following data were obtained: 

1. The deflection curves for the first three bending modes. 

2. The natural bending frequencies corresponding to each of 
these modes. 

3. The center of gravity of the missile. 

4. The mass moment of inertia about the center of gravity. 

5. The total mass of the missile. 


First, the three bending rhodes were plotted as shown in Figures 
2 through 7. The center of instantaneous rotation for each vehicle was 
determined from the relation: 



CG 


mX 


E 


(i) 


where k is the radius of gyration about the vehicle's C. G. and its 
square is equal to the mass moment of inertia about the center of 
gravity, Iq q , divided by the total mass, m. The distance from the 
center of gravity of the vehicle to the swivel point of the engines is 
denoted by 
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The center of pressure, C , is given by the relation: 

P 


C P = 


J^XXx 


dx 


Jfxx 


dx 


(2) 


D 

where X = ■ - and D. . is the diameter of the vehicle as a 

JJ \X/ 

o 

function of x. The diameter at the base of the vehicle is denoted by D . 

o 

The normal lift coefficient, N', was obtained from the relation: 


N» = 7Tq 


•'l 


x 


-f-q D 2 


(3) 


where q is the dynamic pressure. 

The aerodynamic restoring moment coefficient, C^, was calculated 
from the relation? 


7T q° 


c i ■ 


/ 2 X\' ; dJt .. 

I L 


(4) 


C. G. -~ 


The thrust moment coefficient, C^» was obtained from the 
following equation: 

V 


= 


2 I 


( 5 ) 


ca 


where F is the thrust of the swivel engines. Equations 2 through 5 
were obtained from Reference 2. 
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For rigid body control, the equation of motion for the rotation 
of a space vehicle may be written as (c.f. Reference 1) 


vp+ + C 2 j3 = 0 


( 6 ) 


providing translation and sloshing are not present. Assuming no 
wind, a = ; therefore, 

ip + C x f+ C 2 /3 = 0 . (7) 

The simplified control equation can be written as: 


0 = a of + ‘if • 


( 8 ) 


Substituting equation (8) into (7) 

f + C 2 a i V P + [ C i + C 2 a ^^p = 0 • (9) 

It is evident from equation (9) that, for negative, the 
stability of the vehicle is decreased; however, as long as C| + is 
positive, stability is insured providing 0. Since is always 

positive, a^ must be positive also for stability. For the vehicles 
studied in this investigation, all values of are negative. This 
indicates that the normal lift force, N 1 , acts at a point forward to the 
total C. G. of the vehicle and produces a nose-up or destabilizing m 
moment. The fact that the Cp is located forward of the total vehicle 
C. G. is due to the truncated cone and nose cone surfaces which are 
also forward of the C.G. of each vehicle. This aerodynamic instability 
can, however, be eliminated by using engines which can be gimballed 
to produce a stabilizing moment about the vehicle 1 s C. G. The thrust 
moment coefficient C £, then, is a measure of this stabilizing moment 
as defined by equation (5). 
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It can be seen that C 2 /C^ is an indication of the aerodynamic 
stability of the vehicle. The greater C 2 and the smaller C^, the more 
stable the vehicle. To eliminate negative values of the parameter 
C 2 /C 1 ,• the absolute value to will be used as shown in Figure 1. 

Also, the first bending mode frequency divided by the square root of 
the absolute value of is graphed in Figure 1. Since large bending 
frequencies are desirable from the control and sloshing point of view 
and sinc e a small value for Cj is desirable, the vehicle with the larger 
f iA| l c il is the more desirable vehicle, structurally and aerodynamically. 
As a matter of interest, an arbitrary shape parameter (D^ -£-) was 
chosen in an effort to linearize the natural structural bending frequencies 
This parameter is plotted against bending frequencies in Figure 4. The 
weighted diameter of the vehicle, D w , was obtained from equation (10): 


D 

w 


°i ‘2 


D 1 >2 < D 1 + ° 2 > 
+ 7 + 



D 1 
n n 

L 

(10) 


JL 

D 1 

where — is the average diameter of the nose cone, lj is the length 
of the cone, 1^ is the length of the section of vehicle behind the nose 
cone where Dj is constant, etc. Figure 4 enables anyone to determine 
the first three bending frequencies for any vehicle merely be calculating 
the D 4 -j- parameter; however, the weight and internal volume of the 
vehicle must be the same as the Saturn used in this report. 


D 1 *2 (D 1 + D 2 > 
— T"^ + D , *2 + ^2“^ 


1, +. . . D 1 

3 n n 


To determine C^, C^, and N 1 , the maximum dyanmic pressure 
for the C-l trajectory was used, which is 

q MAX = 6. 54 lbs/in^ . 

The thrust from the swivel engines was 


F = 4(188,000) = 752, 000 -lbs . 
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The following table contains the data relevant to this investigation'. 

TABLE I 

Bending and Aerodynamic Stability Data 


Variables* 

Vehicle 

Vehicle 

Vehicle 

Vehicle 

Vehicle 

Vehicle 


0 

i 

2 

3 

4 

5 

f i 

1.237 

2.642 

0.595 

0.308 

0.068 

0.588 

f 2 

2.806 

5.540 

1.498 

0.956 

0.186 

1.533 

f 3 

4.981 

8.409 

2.711 

2.05 

0.361 

2.790 

*IR 

-385.4 

-284.0 

-595.7 

-732.5 

-1149.6 

-500.0 

*MC 

335.7 

151.8 

454.9 

600.9 

54.7 

213.0 

m 

7160 

7101 

7084 

7165 

7150 

7100 

t 

0.210 

0.231 

0.179 

0.126 

0.113 

0.172 

L 

2220 

1730 

2940 

3096 

6798 

3000 

^o* 10 ' 9 

1.86078 

1.23674 

3.88292 

4.44662 

26.66419 

4.5716 

c p 

-655.6 

-700.5 

-401.0 

-481.0 

-3308 

1 

-1393 

N' 

694,100 

694,100 

694,100 

694,100 ! 

147,856 

362,133 

*E 

674.3 

613.3 

920.1 

847.1 

3244.3 

1287 

C 1 

-0.232 

-0.394 

-0.0717 

-0.0752 

-0.0184 

-0.108 

C 2 

0.273 

0.373 

0.178 

0.143 

0.915 

0.211 

y \ c i\ 

1.174 

0.946 

2.486 

1.906 

49.628 

1.955 

f i//W 

2.566 

4.207 

2.220 

1.124 

0.500 

1.789 

i<I R /L 

-0.1736 

-0.1642 

-0.2026 

-0.2366 

j -0.1691 

-0.1667 


195 

225 

172 

162 

116 

177 

A 

rt 

701 

1521 

310 

173 

26 

320 

L 


i 

1 






*Units are in the inch -pound- second system. 

X-j-^ and Cp are measured from the C.G. of the vehicle. 
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The node points for the first bending mode of each missile are 
given in Table II below. The location of each point is measured in 
per cent of total vehicle length aft of the nose. Figures 2 and 3 present 
the vehicle configurations in order of magnitude of their bending 
frequencies and aerodynamic stability decreasing and increasing, 
respectively. 


TABLE II 


NODE LOCATIONS 


VEHICLE 

1st NODE 

2nd NODE 

0 

29. 7% 

79. 7% 

1 

35. 5% 

80. 9% 

2 

30. 6% 

83. 0% 

3 

27. 6% 

83. 2% 

4 

23. 8% 

77. 6% 

5 

28. 3% 

78. 8% 


TYPICAL FIRST BENDING MODE SHAPE 


Yj(X) 


y 

1 . 0 


■ 


^ ^ ^lst NODE 

_ _ - 

NOSE 

— " " "2nd NODE 
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III. CONCLUSIONS 


The graph of / I Cj| shown in Figure 1 indicates that this 
control-stability parameter increases rapidly after the vehicle length 
exceeds 3500 inches. This can be explained by the fact that each of 
the vehicles investigated has the same internal volume; therefor e, its 
diameter decreases with an increase in length. This length increase 
decreases the normal lift force acting on the nose cone in proportion to 
the decrease in the square of base diameter of the vehicle. This causes 
C^ to decrease but' at the same Xp has become larger ; thereby tending 
to increase C^. However, Xp has been increased proportionally to the 
increase in Xp. It should be stated that this is true only for single 
diameter vehicles such as numbers 1, 4, and 5. Since occurs in 
Q 2 and Xp in C^, these two effects cancel and result in C 2 / | C ^ | 
increasing at an exponential rate. 

Although it is desirable for a vehicle to have its first bending 
mode frequency above the control and sloshing frequencies, which can 
be achieved with short vehicles, it is also necessary that any a erody namic 
instabilities be controlled by swivel engines. The graph of f ]/ v/7^ll 
(Figure 1) illustrates the variation of this parameter with vehicle length 
and allows one to choose, when compared with the C^/ )C]J parameter, 
some compromise design between maximum bending frequency-minimum 
stability, and maximum stability-minimum bending frequency. Vehicles 
1 and 4 are the extreme cases investigated here. For conical, single- 
diameter cylindrical, and dual -diameter cylindrical vehicles with the 
same approximate lengths, e. g. , vehicles 2, 5, and 3, respectively , the 
conical vehicle exhibits more favorable bending and control characteristics 
than the other two, while the single-diameter cylindrical is better than 
the dual-diameter cylindrical vehicle. 

Figure 4 illustrates the variation of the first three bending modes 
of each vehicle with the shape parameter, L. A linear relationship 
exists for the first mode for all values of shap‘d parameter while the 
second mode is linear only for values of D ^ L greater than 400. The 
third mode is nearly linear for D ^ gr eater-'than 500. The bending 
frequencies obtained for each vehicTe are lower than those of an actual 
vehicle designed to the length and diameters specified in this report. 

This is due mainly to the way that the El and mass distributions were 
obtained (see Description). The bending frequency of vehicle number 
zero is about 70 per cent of the actual frequency of SA-1 at lift-off. 
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It has been found that, when the propellant sloshing mass is 
located between the center of gravity and the center of instantaneous 
rotation of a vehicle, its stability is endangered; therefore, it is 
desirable to have these two positions located as close to each other 
as possible. Generally, the sloshing mass will be located aft of the 
vehicle center of gravity for first stage flight time and the center of 
instantaneous rotation forward of the C. G. From Table I it is seen 
that the more stable vehicles from this standpoint are, in descending 
order: 1, 5, 4, 0, 2, and 3. This is evidenced by the magnitude of 
Xir/L. In general, the smaller this value, the less the probability 
that an instability will occur due to sloshing, provided that the vehicle 
has only one sloshing mass, i. e.^that sloshing masses of all other 
tanks are small in comparison to this one. It is also necessary that 
the sloshing masses be equal in each vehicle for the above statements 
to be true, which is not the case in the actual vehicle. For the small 
diameter vehicles, the sloshing mass is small, while for vehicles with 
large diameters the sloshing mass is dangerously large. 




FIGURE I. FREQUENCY AND AERODYNAMIC STABILITY PARAMETER 
CURVES FOR THE VARIOUS CONFIGURATIONS 
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FIGURE 2. VEHICLE CONFIGURATIONS IN ORDER OF THEIR 
BENDING FREQUENCIES (DECREASING) 









FIGURE 4. BENDING FREQUENCY VS SHAPE PARARMETER 
CURVES FOR THREE MODES 
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APPENDIX A 


The first three bending mode deflection curves for each 
vehicle are given in this appendix. 



FIGURE 5. BENDING MODE DEFLECTION CURVES 
FOR CONFIGURATION #0 
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FIGURE 6. BENDING MODE DEFLECTION CURVES 
FOR CONFIGURATION #1 





FIGURE 7. BENDING MODE DEFLECTION CURVES 
FOR CONFIGURATION # 2 







FIGURE 8. BENDING MODE DEFLECTION CURVES 










22 



FIGURE 10. BENDING MODE DEFLECTION CURVES 





APPENDIX B 


The mass distributions both empty and full for each 
configuration are contained in this appendix. 



MASS DISTRIBUTION FOR CONFIGURATION *0 
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FIGURE 12. MASS DISTRIBUTION FOR CONFIGURATION # I 
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FIGURE 13. MASS DISTRIBUTION FOR CONFIGURATION #2 







14. MASS DISTRIBUTIONS FOR CONFIGURATION #3 







EMPTY FULL 
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FIGURE 15. MASS DISTRIBUTION FOR CONFIGURATION #4 






empty Fy LL 
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